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The study of Bhabha scattering at e+e — colliders probes the running of the electromagnetic coupling. After 
early measurements by the VENUS collaboration at TRISTAN and the by L3 collaboration at LEP, two recent 
analyses have been performed by the OPAL and L3 collaborations. The OPAL collaboration studied high- 
statistics low-angle Bhabha scattering at LEP, achieving a precise determination of the running of a in the region 
1.8 GeV 2 < — Q 2 < 6.1 GeV 2 . The L3 collaboration investigated high-energy large-angle Bhabha scattering to 



first probe the region 1800 GeV 2 < -Q 2 < 21600 GeV 2 . 
of their agreement with QED predictions is given. 



All measurements are described and a global overview 
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1. THE RUNNING OF a 

A fundamental idea in quantum field theory is that 
coupling constants which describe the strength of a 
given process are not actually constant but rather de- 
pend on the energy scale at which the process occurs; 
in other words, the constants run. In QED this run- 
ning implies the increase of the electromagnetic cou- 
pling, a, with the squared four-momentum transfer, 
Q 2 . Figure 1 sketches the origin of this phenomenon in 
the case of Bhabha scattering, e + e~ — > e + e~: larger 
momentum transfers probe virtual-loop corrections to 
the photon propagator. 

The running of a is parametrised as [1]: 



a(Q 2 ) = 



Of) 



1 - Aa(Q 2 



(1) 



where olq is the fine-structure constant. This is mea- 
sured with high accuracy in solid-state processes and 
via the study of the anomalous magnetic moment of 
the electron to be [2]: 

l/a = 137.03599911 ± 0.00000046 

The term Aa(Q 2 ) describes the running of a in Equa- 
tion 1. It receives contributions from three sources: 
lepton loops, top-quark loops and loops involving the 




Figure 1: i-channel Feynman diagrams contributing to 
Bhabha scattering and the phenomenon of 
vacuum-polarisation. The sum of all diagrams including 
zero, one, two or more vacuum-polarisation insertions is 
denoted by the diagram to the left with the double-wavy 
photon propagator, with an electromagnetic coupling 
a(Q 2 ). 



five lighter quark flavours. At the scale of Z-boson 
mass, mz, the first two contributions to Aa(Q 2 ) are 
precisely calculated as: 



0.03150 



Aa top (m|) = -0.00007 ±0.00001, 



(2) 



(3) 



as discussed in References [3] and [4], respectively. 
The contribution from the five lighter quarks, Aa^, 
is difficult to calculate due to non-perturbative QCD 
effects and is estimated using dispersion-integral tech- 
niques as [5]: 



Aa hL( m z) = 0.02758 ± 0.00035. 



(4) 



The electromagnetic coupling at the scale of the Z- 
boson mass is therefore [5]: 



l/a(m|) = 128.940 ±0.048. 



(5) 



The 



study of fermion-pair production at LEP, 
c + e~ — > Z — > ff , allows a precise measurement of the 
Z-boson couplings. These couplings are sensitive to 
electroweak vacuum-polarisation effects and therefore 
also to the running of a. The quantity most sensitive 
to the running of a is the effective vector coupling 
constant of Z bosons and electrons. It is measured 
as g Vc = -0.03816 ± 0.00047 [6]. In the absence of 
the running of a its value would be gy c — —0.076. A 
global fit to several Standard Model observables al- 
lows to constrain the value of a at the scale of the 
Z-boson mass as [6] : 



l/a flt {ml) = 128.937 ± 0.047, 



(G) 



which is in remarkable agreement with the expected 
value of Equation 5. However, this determination of 
a(m|) is based on the assumption of the full elec- 
troweak theory and a determination which relies on 
less hypotheses is desirable. The study of Bhabha 
scattering, discussed in the following, allows a more 
direct insight in the running of a. 
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2. BHABHA SCATTERING AND THE 
RUNNING OF a 

The study of Bhabha scattering at e + e~ colliders 
offers a unique window on the vacuum-polarisation in- 
sertions depicted in Figure f. Its measurement gives 
access to the running of a in the space-like region, 
Q 2 < 0. The experimental conditions at e + e~ col- 
liders are extremely favourable as they allow a precise 
determination of the four-momentum transfer through 
its dependence on the squared centre-of-mass energy, 
s, and on the scattering angle, 9: 

Q 2 = t~ -s(l-cos0)/2. 

Indeed, at e + e~ colliders, s is accurately known and 
detectors are conceived for a precise measurement of 
9. In the following, Bhabha scattering is studied in 
two angular regions: low-angle, 9 <~ 1.5° — 3.5°, and 
large-angle, 9 ~ 20° — 90°. These angular regions 
correspond to two energy regimes: low-Q 2 and high- 
Q 2 , respectively. 

The relation between the differential cross section of 
Bhabha scattering and the electromagnetic coupling 
is: 



dcr dcr° / aft) , , 



where the tree-level cross section is: 



z, 



da^ 
dt 



Anal 



(7) 



(8) 



The s-channel contributions, <5 7 and 5%, are much 
smaller than the radiative corrections, e. These three 
quantities are known with a good accuracy compared 
to the experimental precision of the study of low- 
and large-angle Bhabha scattering. Naively, one could 
imagine to measure the cross section for Bhabha scat- 
tering and then insert the result in the left-hand side 
of Equation 7 in order to extract a value for a(t). 
This argument is unfortunately flawed: the measure- 
ment of the cross section requires knowledge of the 
integrated luminosity. At LEP this is estimated by 
counting events from low-angle Bhabha scattering and 
assuming the cross section of this process to be known, 
what can only happen if a(t) is known to start with: 
a catch-22 [7]! 

As discussed above, the information on a(t) con- 
tained in the scale of the differential cross section 
da/dt cannot be directly unlocked. On the other 
hand, the shape of this differential cross section con- 
tains information on the running of a(t) over the Q 2 
range in which the cross section is measured. QED 
predicts that a(t) increases with t: if the running is 
slower than the one expected in QED, the differen- 
tial cross section will be steeper than expected; if the 
running is faster than expected in QED, the differen- 
tial cross section will be flatter than expected. The 



strategy to access information on the running of a is 
therefore to analyse the shape of the differential cross 
section over a Q 2 range and compare it with the ex- 
pectation from QED. Five such analyses exist and are 
presented below: we first discuss three early measure- 
ments of the running of a at TRISTAN and LEP, and 
then describe two recent measurements performed by 
the OPAL and L3 collaborations by investigating low- 
angle and large-angle Bhabha scattering at LEP. 



3. EARLY MEASUREMENTS AT TRISTAN 
AND LEP 



LEP experiments were equipped with low-angle 
calorimeters, called luminosity monitors, engineered 
to detect events from low-angle Bhabha scattering in 
order to measure the integrated luminosity. The event 
counts in different polar-angle regions of the luminos- 
ity monitors depend on the shape of the differential 
cross section of Bhabha scattering and, through Equa- 
tion 7, give access to the running of a(t). 

The first study of the running of a using low-angle 
Bhabha scattering was performed by the L3 collabo- 
ration. A total of about 7 million events collected at 
\fs = mz in the luminosity monitor, covering a polar 
angle 1.8° < 9 < 3.1°, were used to probe the running 
of a in the region 2.1 GeV 2 < -Q 2 < 6.2 GeV 2 [8]. 
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Figure 2: Ratio of measured and expected events in four 
different angular regions of the L3 luminosity monitor, 
corresponding to four values of — Q 2 , for a subset of data 
collected at the Z-boson pole. The solid circles 
correspond to the QED prediction and the open circles to 
the hypothesis of no running of the electromagnetic 
coupling. The data favours the QED predictions and 
excludes the no-running hypothesis. From Reference [8] . 



Figure 2 presents the ratio between the observed 
number of events and the QED prediction, as well 
as the hypothesis of a constant value of a over the 
entire Q 2 range. This ratio of data and predictions is 
fitted to determine the running of a between the two 
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Figure 3: Ratio of the differential cross section for 
Bhabha scattering measured by the VENUS 
collaboration and the predictions for a fixed value of a. 
The solid line shows the QED prediction for a running a 
and the dashed line the prediction if the running is only 
due to lepton corrections. The data favours the QED 
predictions. From Reference [9]. 

extremes of the Q 2 range as: 

a- 1 (-2.1 GeV 2 ) -a _1 (-6.2 GcV 2 ) = 0.78 ± 0.26, 

(9) 

where the uncertainty combines statistical and sys- 
tematic uncertainties. This result is in good agree- 
ment with the QED prediction of 0.56 and first estab- 
lished the running of a in this energy range, excluding 
the no-running hypothesis at 4.5 standard deviations. 

The running of a in large-angle Bhabha scattering 
was first investigated by the VENUS Collaboration 
at TRISTAN. They studied events collected at y/s = 
57.8 GeV in the angular range 47° < 6 < 90° to cover 
the Q 2 range 100 GeV 2 < -Q 2 < 3337 GeV 2 [9]. 

Figure 3 shows the measured differential cross sec- 
tion as a function of the cosine of the scattering angle, 
normalised to the QED predictions with a fixed value 
of a. The effects of the running are evident and the 
data favour a scenario in which both the hadronic and 
leptonic contributions are present in the running, as 
opposed to a scenario with just lepton loops. The 
results can be expressed as the change over the Q 2 
range of the term Aa describing the running of a in 
Equation 1: 

Aa(-3337 GcV 2 ) - Aa(-100 GcV 2 ) = 

0.0186 ±0.0059, (10) 

where the uncertainty includes statistical and system- 
atic effects. 

A similar study at higher centre-of-mass energies 
and lower polar angles was performed by the L3 Col- 
laboration. Bhabha-scattcring events are selected at 
y/s = 189 GeV for scattering angles 20° < 6 < 
36°, probing the intermediate-Q 2 range 12.25 GcV 2 < 
-Q 2 < 3434 GeV 2 [8]. The measured differential 
cross section as a function of the scattering angle is 



compared to the predictions as a function of a, allow- 
ing to measure its running over the Q 2 range as: 

a' 1 (-12.25 GeV 2 ) - a~ 1 (-3434 GeV 2 ) = 3.8 ± 1.3, 

(11) 

where the uncertainty comprises statistical and sys- 
tematic sources. This value is in good agreement with 
the QED prediction of 4.07. 

4. PRECISION MEASUREMENT AT LOW 

Q 2 

The OPAL collaboration analysed low-angle 
Bhabha scattering events to study the running of 
a at low Q 2 [10]. Data collected at yfs = mz 
with the luminosity monitor, consisting of layers 
of tungsten absorber and 32-pad silicon detectors, 
were used. This detector covered a polar region 
1.4° < 6 < 3.3°, corresponding to a momentum 
transfer 1.8 GcV 2 < -Q 2 < 6.1 GcV 2 . 

The analysis selects 10 million high-energy back- 
to-back cluster pairs and studies their t spectrum to 
extract information on a(t). Data are divided in five 
t bins and are compared with the BHLUMI Monte 
Carlo [11]. The results of the study are quantified 
by studying the ratio of the t spectrum observed in 
data and the Monte Carlo spectrum for the hypothesis 
a(t) = cto, as shown in Figure 4. 
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Figure 4: The ratio of event counts five radial regions of 
the OPAL luminosity monitor, corresponding to five Q 2 
values, for data and a Monte Carlo simulation with 
a(t) = ao- The data favour the QED prediction for the 
running of Aa(t) and exclude scenarios with no running 
or a running due only to electron or lepton virtual loops. 
Only statistical uncertainties are shown. From 
Reference [10]. 

The data in Figure 4 are fitted with the function 
a + bln(t/t ), where to = —3.3 GeV 2 is the mean value 
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of t in the data sample. The parameter b is related to 
the running of a by: 



Aa(t 2 ) - Aa(ii) 



b. h 
- in — . 

2 t 1 



(12) 



Four hypotheses are tested by the OPAL analysis: no 
running of a; a running only due to electron virtual- 
loops; a running only due to lepton virtual-loops; the 
QED hypothesis of a running induced by virtual loops 
of lcptons and quarks. The first three hypotheses are 
found to be completely excluded by data. 

The observed running of a is derived from the fitted 
value of b and Equation 12 to be: 



Aa(-6.1 GeV 2 ) - Aa(-1.8 GcV 2 ) = 
(440 ± 58 ± 43 ± 30) x 10" 5 , 



(13) 



where the first uncertainty is statistical, the sec- 
ond systematic and the third theoretical. The main 
sources of systematic uncertainties are the simulation 
of the detector material and the reconstruction of 
the radial coordinates of the energy clusters in the 
calorimeter. This results is in excellent agreement 
with the QED prediction of 460 x 10" 5 . 

This measurement proves the running of a at low 
Q 2 with a significance of 5.6a. In addition, it estab- 
lishes for the first time the hadronic contribution to 
the running with a significance of 3. Oct as: 

Aa had (-6.07 GcV 2 ) - Aa had (-l.Sl GcV 2 ) = 

(237 ± 58 ± 43 ± 30) x 10~ 5 . (14) 



5. FIRST MEASUREMENT AT LARGE Q 2 

The L3 collaboration probed the running of a at 
large Q 2 by studying large-angle Bhabha scattering 
in the angular region 26° < 9 < 90° [12]. Data col- 
lected at eight centre-of-mass energies in the range 
= 189 — 209 GeV were considered, corresponding 
to a four-momentum transfer 1800 GeV 2 < — Q 2 < 
21600 GeV 2 . About 40 000 events are selected as 
back-to-back clusters in the high-resolution BGO elec- 
tromagnetic calorimeter with matched tracks. These 
events are used to measure the differential cross sec- 
tion of Bhabha scattering in ten different angular 
ranges, for each of the eight energy points, for a total 
of 80 measurements. The measured cross section are 
compared to the predictions of the BHWIDE Monte 
Carlo [13] to extract information on the running of a. 
As an example, Figure 5 compares the data and the 
predictions for the ten centre-of-mass-averaged cross 
sections. 

The results of the study are expressed by inserting 
an additional parameter, C, in the description of the 
running of a of Equation 1 : 



a(Q 2 ) 



a 



^ 1.2 

<x> 
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o 
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Figure 5: The ratio of the centre-of-mass-averaged 
differential cross section of for large-angle Bhabha 
scattering measured by L3 and the corresponding 
theoretical predictions. The inner error bars denote the 
statistical uncertainties, the outer the combination of 
statistical and systematic uncertainties. The data 
exclude the scenario a(t) = ao, represented by the dotted 
line. From Reference [12]. 



A fit to the 80 data points, which takes into account 
the dependence of the measured integrated luminosity 
on C shows that the value C = 0, corresponding to 
absence of running, is completely excluded. The data 
are in excellent agreement with the running predicted 
in QED, corresponding to C = 1, and the fit yields: 



C = 1.05 ±0.07 ±0.14, 



(16) 



1 - CAa(Q 2 ) ' 



(15) 



where the first uncertainty is statistical and the second 
systematic. The systematic uncertainty is dominated 
by theoretical uncertainties on the prediction of the 
Bhabha scattering differential cross section both at 
large angles and in the luminosity monitor. Some ad- 
ditional systematic contributions also come from the 
modelling of the detector response. 



6. COMBINED RESULTS 

Figure 6 summarises the LEP results on the run- 
ning of the electromagnetic coupling. The L3 mea- 
surement at 1800 GeV 2 < -Q 2 < 21600 GcV 2 [12] is 
represented as a band, obtained by inserting the mea- 
sured value of C from Equation 16 into Equation 15 
and assuming the QED description of Aa(Q 2 ) given in 
Reference [5]. The two low-Q 2 measurements by the 
OPAL and L3 collaborations at 1.8 GcV 2 < -Q 2 < 

6.1 GeV 2 (Equation 13) and 2.1 GcV 2 < -Q 2 < 

6.2 GeV 2 (Equation 9), respectively, as well as the 
intcrmcdiatc-Q 2 measurement of L3 at 12.25 GcV 2 < 
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Figure 6: Summary of LEP results on the measurement 
of the running of the electromagnetic coupling. The band 
represents the L3 measurement at high Q 2 . The full 
symbols represent the OPAL and the L3 measurements 
at low and intermediate Q 2 . The open symbols are the 
reference values to which the measurement are anchored, 
as discussed in the last section of the text. The solid line 
shows the QED predictions of Reference [5] . 



-Q 2 < 3434 GeV 2 (Equation 11) are represented with 
two symbols each. The empty symbols represent the 
values of a^ 1 (Q 2 ) at the lower end of each Q 2 range. 
They are fixed by using Equation 1 and the QED de- 
scription of Aa(Q 2 ) of Reference [5]. The full symbols 
represent the values of a~ 1 (Q 2 ) at the higher end of 
each Q 2 range extracted from the fixed values and 
from the measurements in Equations 13, 9 and 11. 
All measurements are in excellent agreement with the 
QED predictions of Reference [5] , presented as a solid 
line. 

Figure 7 presents a combination of all LEP results. 
This combination is obtained in several steps. First, 
the L3 measurement at 2.1 GeV 2 < -Q 2 < 6.2 GeV 2 
and the OPAL measurement at 1.8 GeV 2 < — Q 2 < 
6.1 GeV 2 , from Equations 13 and 9 respectively, are 
combined into a single measurement. In order to per- 
form this combination the L3 measurement is trans- 
ported to the Q 2 range of the OPAL measurement. 
The combined result is: 

a(-6.1 GeV 2 ) - a(-1.8 GeV 2 ) = (363 ± 52) x 1CT 7 , 

(17) 

where the uncertainty combines statistical and sys- 
tematic effects. This combined result and the L3 
measurement at intermediate Q 2 are plotted in Fig- 
ure 7 as two single points at Q 2 — —6.1 GeV 2 and 
Q 2 = -3434 GeV 2 , respectively. The values of a(Q 2 ) 
at these two points are extracted by anchoring the 
value of a(Q 2 ) at the lower end of each Q 2 range 
by using the L3 measurement of C at 1800 GeV 2 < 
-Q 2 < 21600 GeV 2 and assuming it also describes 



the running of a for lower values of Q 2 : the values of 
a(-1.8 GeV 2 ) and a(-12.25 GeV 2 ) are fixed by us- 
ing the measured value of C from Equation 16, the 
evolution expected from Equation 15 and the QED 
description of Aa(Q 2 ) of Reference [5]. The value of 
a(— 6.1 GeV 2 ) is finally extracted by using this fixed 
value of a(— 1.8 GeV 2 ) and Equation 17, with an addi- 
tional uncertainty which follows from the 14% uncer- 
tainty on C. A similar procedure is followed to extract 
the value of a(— 3434 GeV 2 ) from Equation 11. 
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Figure 7: Combined LEP results on the measurement of 
the running of the electromagnetic coupling in three 
different Q 2 regimes, compared with the QED 
predictions of Reference [5] . The treatment of data is 
discussed in the last section of the text. 

In conclusion, the LEP experiments have estab- 
lished the the evolution of the electromagnetic cou- 
pling with the squared four-momentum transfer in a 
new energy domain. These measurements, combined 
in Figure 7, span three orders of magnitude in Q 2 and 
confirm the QED predictions for the running of a. 
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